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Recent results on the spectral properties of the Hermitian Wilson-Dirac operator are presented. 



1. Introduction 

Continuum gauge field theory works under the 
assumption that all fields are smooth functions 
of space time. This assumption is certainly a 
valid one for quantum gauge field theories that 
respect gauge invariance: One should always be 
able to fix a gauge so that the gauge fields are 
smooth functions of space time since the action 
that contains derivatives in gauge fields will not 
allow it otherwise. The space of smooth gauge 
fields typically has an infinite number of discon- 
nected pieces with the number of pieces being in 
one to one correspondence with the set of inte- 
gers 0. Every gauge field in each piece can be 
smoothly interpolated to another gauge field in 
the same piece but there is no smooth interpola- 
tion between gauge fields in different pieces. This 
is the case for U(l) gauge fields in two dimensions 
and SU(N) gauge fields in four dimensions. 

In lattice gauge theory, gauge fields are repre- 
sented by link variables U fJi (x) that are elements 
of the gauge group. Continuum derivatives are 
replaced by finite differences and the concept of 
smoothness of gauge fields does not apply. Any 
lattice gauge field configuration, U^x) — e %A ^ x ^ 
can be deformed to the trivial gauge field config- 
uration by the interpolation J7 M (x;r) = e lTA ^^ 
with U^x]!) — U f _ l (x) and [/ M (a;;0) = 1. Since 
smoothness does not hold on the lattice away 
from the continuum limit, the space of gauge 
fields on the lattice forms a simply connected 
space. Separation of the gauge field space into an 
infinite number of disconnected pieces can only 
be realized in the continuum limit. 

In this talk, I will address the following basic 



question: Do we see a separation of lattice gauge 
fields configurations into topological classes as we 
approach the continuum limit? To answer this 
question, I will assume that an ensemble of lat- 
tice gauge field configurations are given to me. 
This could be a pure gauge field ensemble or an 
ensemble that includes fermion dynamics using ei- 
ther staggered fermions or Wilson fermions. My 
aim is to answer this question without altering 
the lattice gauge field configuration. I will use a 
Wilson-Dirac fermion to probe the lattice gauge 
field configuration. My motivation is the over- 
lap formalism Q for chiral gauge theories. Topo- 
logical aspects of the background gauge field are 
properly realized by the chiral fermions in this 
formalism and therefore it provides a good frame- 
work to answer the above question. The hermi- 
tian Wilson Dirac operator enters the construc- 
tion of lattice chiral fermions in the overlap for- 
malism and topological properties of the gauge 
fields are studied by looking at the spectral flow 
of the hermitian Wilson Dirac operator as a func- 
tion of the fermion mass. 

The talk is organized as follows. I begin by ex- 
plaining the connection between the spectral flow 
of the hermitian Wilson Dirac operator and the 
topological content of the background gauge field. 
I then present possible scenarios for the qualita- 
tive nature of the spectrum on the lattice. I will 
then present some properties of the spectral flow, 
and I will illustrate these properties by studying 
the spectral flow in a single instanton background. 
I will then focus on the spectral properties on lat- 
tice gauge field ensembles and their behavior as 
the continuum limit is approached. I will only 
present an overview of the results. Detailed re- 
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suits can be found in the talk presented by Ed- 
wards 101. 



2. Spectral flow, topology and condensates 

The massless Dirac operator in the continuum 
anticommutes with 75. Therefore, the non-zero 
imaginary eigenvalues of the massless Dirac op- 
erator come in pairs, ±zA with ip and "f^ip being 
the two eigenvectors. The zero eigenvalues of the 
massless Dirac operator are also eigenvalues of 75 . 
These chiral zero modes are a consequence of the 
topology of the background gauge field. It is use- 
ful to consider the spectral flow of the Hermitian 
Dirac operator: 

H(m) =7&(7/*A. + "*) (1) 

The non-zero eigenvalues of the massless Dirac 
operator combine in pairs to give the following 
eigenvalue equation: 

H(m)x± = ± VA 2 + m 2 x± (2) 

X± are linear combinations of ip and 75^ and 
these modes never cross the x-axis in the spec- 
tral flow of H(m) as a function of m. The zero 
eigenvalues, 7 M -D M 0± = with 75</>± = ±4>± re- 
sult in 

H(m)4>± = ±m<f>± (3) 

These modes associated with topology result in 
flow lines that cross the x-axis. A positive slope 
corresponds to positive chirality and vice-versa. 
The net number of lines crossing zero (the dif- 
ference of positive and negative crossings) is the 
topology of the background gauge field. Global 
topology of gauge fields cause exact zero eigen- 
values at m = 0. In addition, one can have a 
non-zero spectral density at zero. At infinite vol- 
ume in the continuum, the spectrum is continuous 
and p(A; m)d\ is the number of eigenvalues in the 
infinitesimal region dX. The spectral gap A 9 (m) 
defined as the lowest eigenvalue at to is equal to 
m|. The spectral density at zero, p(0;m) can be 
non-zero only at m = indicating chiral symme- 
try breaking in a theory like QCD. 

To study the possible emergence of the above 
continuum picture as a continuum limit of a lat- 
tice gauge theory picture, we need to have a lat- 
tice realization of H(m). It is important to note 
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Figure 1. Continuum picture of the spectral gap 
and the spectral density at zero. 



that we are interested in the spectral flow of a sin- 
gle Dirac fermion. With this in mind, we choose 
the hermitian Wilson-Dirac operator obtained by 
multiplying the standard Wilson-Dirac operator 
by 75: 



Hi (to) = 



B(U) +r 
Ct(tf) 



C(U) 
-B(U) - 



(4) 



if to < 



Let 



be a normalized eigen- 



C is the naive lattice first derivative term and 
B is the Wilson term. We are interested in the 
spectral flow of H^(m) as a function of m. We 
note that to = 0, —2, —4, —6, —8 are the points 
where the free fermions become massless with de- 
generacies 1,4,6,4,1 respectively. Next we ob- 
serve that Hz, (to) can have a zero eigenvalue only 
u 
v 

vector of Hi(m) with zero eigenvalue. Then 
mu + Bit + Gv — and C^u — Bv — mv = 0. 
This implies that u^Bu + v'Bv = —to. Since B 
is a positive definite operator, the above equation 
can have a solution only if to < 0. 

We focus on the range — 2 < m < and pro- 
pose the following scenarios for the spectral gap 
and the spectral density at zero on the lattice and 
their approach to the continuum limit. Six differ- 
ent scenarios are possible as shown in Fig. ^. 

• On the lattice we have (a) and (i) with 
to c — > in the continuum limit. p(0; 0) ap- 
proaches the continuum limit with proper 
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Figure 2. Possible scenarios of the spectral gap 
and the spectral density at zero on the lattice. 



scaling taken into account. 

• On the lattice we have (b) and (ii) with (a) 
and (i) being the continuum limit. In this 
case, mi — > and m% — ► 0. In the limit we 
also get p(0; 0). 

• On the lattice we have (c) and (ii) with (a) 
and (i) being the continuum limit. In this 
case, mi — > 0. The gap opens up for m < 
and again we get p(0; 0). 

• On the lattice we have (c) and (ii) with (c) 
and (i) being the continuum limit. In this 
case, mi — > 0. But the gap does not open 
up for m < 0. 

• On the lattice we have (c) and (iii) with (a) 
and (i) being the continuum limit. In this 
case, mi — > 0. The gap opens up for m < 
and again we get p(0; 0). 

• On the lattice we have (c) and (iii) with (c) 
and (i) being the continuum limit. Here also 



mi — » and p(0, m) = if m < 0. But the 
gap does not open up for m < 0. 

I will show that numerical studies of the spectral 
flow on various ensembles favor the last scenario. 
This seems to be the case for any lattice gauge en- 
semble whether the ensemble was generated with 
fermion dynamics or not. One should keep in 
mind that the Wilson fermion used to the study 
the spectral flow is only a probe and does not en- 
ter the dynamics that generated the gauge field 
ensemble. 

3. Properties of the lattice spectral flow 

I now present a topological argument which will 
show that zero eigenvalues of Hi(m) can occur 
anywhere in the region — 8 < m < || . The spec- 
trum of Hi(m) and — Hl(— 8— m) are identical for 
an arbitrary gauge field background. Since zero 
eigenvalues can occur only for m < in Hi(m), 
it follows that zero eigenvalues can occur only in 
the region — 8 < m < 0. It also follows that every 
level crossing zero from above in the spectral flow 
of (m) should be accompanied by a level cross- 
ing zero from below. In a smooth instanton back- 
ground a level crossing zero from above at mi is 
accompanied by another level crossing zero from 
below at ni2 < m\. The second crossing is due to 
one of the four doubler modes. Both mi and ni2 
will be functions of the size of the instanton p in 
lattice units. For p >> a, mi w and m.2 ~ —2. 
As p decreases, mi moves farther away from zero 
and m,2 moves away from —2 and closer to mi. 
This motion as a function of p is smooth and for 
some value of p, mi = mi- Spectral flow changes 
smoothly as the configuration is slowly changed. 
As we move in configuration space the topolog- 
ical charge of a configuration changes. Tracing 
the spectral flow as a function of configurations 
shows that zero eigenvalues of Hi(m) can occur 
anywhere in the region —8 < m < 0. 

Fig. H illustrates the various properties of the 
spectral flow of the Hermitian Wilson-Dirac op- 
erator. The background gauge field is a smooth 
instanton with p = 2.0 on a 8 4 lattice. The first 
crossing on the right of the figure from above to 
below is the one associated with the particle feel- 
ing the global topology of the gauge field. Then 
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there are four levels crossing zero out of which 
three are degenerate. The degeneracy is a conse- 
quence of the smoothness of the lattice instan- 
ton. The four crossings are attributed to the 
four doublers feeling the global topology. This 
is followed by six levels crossing zero, four levels 
crossing zero and finally one level crossing zero. 
Clearly the spectrum of Hz, (to) and — Ht(8 — m) 
are the same. There is no net level crossing in 
the range — 8 < to < 0. This is because the chi- 
ralities of the sixteen particles are paired. Fig. || 
shows the location of the first level crossing zero 
as one varies the size of the instanton in lattice 
units. As the size of the instanton is decreased 
the crossing point moves farther away from zero. 
p < 1 instantons are placed on the lattice by cen- 
tering the instanton on a lattice site. This figure 
shows that one can move the first crossing all the 
way up to to = —2. 



Figure 3. Spectral flow of a smooth SU(2) instan- 
ton with p = 2.0 on an 8 4 lattice 



Figure 4. toi versus p for a smooth instanton 



4. Spectral density at zero 

In the previous section, I argued that Hi(m) 
can have zero crossings anywhere in the region 
—2 < to < 0. Therefore the spectral gap is zero 
in this region on the lattice. This has direct im- 
plications for how the spectral density at zero 
behaves on the lattice. A careful study of the 
spectral density at zero has been performed on a 
variety of SU(3) pure gauge ensembles ||. This 
is done by computing the low lying eigenvalues 
of Hl (to) using the Ritz functional Q|. The low 
lying eigenvalues over the whole ensemble is then 
used to obtain the integral of the spectral density 
function, namely f Q p(X';m)d\ / . This is then fit- 
ted to a polynomial and p(0; to) is obtained as 
the coefficient of the linear term in the fit. All 
ensembles show a peak in p(0; m) at some value 
of to. There is a sharp rise to the peak from the 
right and a gradual fall from the peak on the left. 
There is a gradual rise again to a second peak at 
the location of the first set of doublers. The peak 
itself gets sharper as one goes toward the contin- 
uum limit. It also moves to the right. p(0;m) is 
non-zero for — 2 < to < toi in the infinite vol- 
ume limit at any finite value of the lattice gauge 
coupling, toi goes to zero as the lattice coupling 
approaches the continuum. p(0; to) approaches 
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Figure 5. p(0;m) as a function of to for various 
SU(3) pure gauge ensembles. 



the infinite lattice volume limit from below as ex- 
pected. We are fairly confident that we have the 
infinite volume limit estimate for p(0; m) at all 
the lattice spacings plotted in Fig. |[ In Fig. |[ 
I focus on the behavior of p(0; m) at a fixed to 
as one approaches the continuum limit. The first 
figure indicates that p(0; m) goes to zero expo- 
nentially in the inverse lattice spacing. This is 
given some credence by plotting the same figure 
in a logarithmic scale. The function seems to fa- 
vor a fit of the form be~ c ^ v". The power of a in 
the exponent is a consequence of an emperical fit 
but the data presents substantial evidence for the 
following: p(0; m) is non-zero for all finite lattice 
spacings. The approach to zero at zero lattice 
spacing is faster than any power of the lattice 
spacing. This shows that the last scenario pre- 
sented in the previous section is favored by the 
numerical results on the lattice. 

5. Discussion 

A probe of lattice gauge fields using Wilson 
fermions has revealed that fields are not contin- 
uum like on the lattice at gauge couplings that 
are typically considered to be weak. If they were 
continuum like, we should have seen evidence that 
p(0, to) is non-zero at a single value of to or a re- 
gion in to that is of the order of the lattice spac- 
ing. Further we should have seen a symmetry 
in the spectrum at values of to on cither side of 
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Figure 6. The approach of p(0; —1.7) to the con- 
tinuum limit 



the point (or region) where p(0, to) is non-zero. 
Instead, we found that p(0, to) is non-zero in a 
region —2 < m < mi- In the continuum limit, 
there is evidence that toi goes to zero and that 
p(0, to) goes to zero away from m = 0. But the 
spectral distribution docs not show evidence for 
a symmetry as m — > —to. 

In addition to studying p(0, to), one should also 
look at the density of levels crossing zero in an 
infinitesimal range dm centered at to. In the con- 
tinuum we expect level crossings only at to = 
but we find a finite density of level crossings zero 
where ever p(0;m) is non-zero ||. This implies 
that the topological charge of a lattice gauge field 
configuration defined as the net level crossings 
in Hi(TO) in the range [to,0] will depend on to. 
The topology of a single lattice gauge field con- 
figuration is not interesting in a field theoretic 
sense. One has to obtain an ensemble average of 
the topological susceptibility and study its depen- 
dence on to. This has been studied on a variety 
of ensembles |||| and the results show that the 
topological susceptibility is essentially indepen- 
dent of to in the region to the left of the peak in 
p(0;m). This indicates that the levels that cross 
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zero to the left of the peak are not physically 
relevant consistent with the result that p(0; to) 
goes to zero in the continuum limit. The numbers 
for the topological susceptibility obtained in this 
manner are consistent with the results obtained 
by field theoretic methods ||. The modes that 
cross zero to the left of the peak in p(0; to) corre- 
spond to small modes ||[| and this is consistent 
with Fig. |. 

A remark on the approach of p(0, to) to the 
thermodynamic limit at a fixed lattice gauge cou- 
pling is in order. For definiteness, let us consider 
a pure gauge ensemble. If the extent of the lat- 
tice is smaller than the extent needed for a finite 
temperature phase transition then p(0, to) will be 
essentially zero. One will see a sharp increase in 
p(0, to) at lattice volumes close to the phase tran- 
sition. A table of f3 c (N T ) can be found in Q for 
SU(3) and SU(2). For example, one needs to have 
a lattice volume bigger than 8 4 if one is using an 
SU(3) gauge coupling of (3 = 6.0 and one needs to 
work with a lattice larger than 14 4 at (3 — 6.38. 
Only then will one see a value of p(0, to) close to 
the thermodynamic limit. 

Studies of the spectrum of the Wilson-Dirac 
operator is not in any sense a new field of re- 
search. Smit and Vink jl0| studied the complex 
spectrum of the Wilson-Dirac operator to under- 
stand topology. These initial studies led to the 
understanding of dislocations in lattice field con- 
figurations that affect the field theoretic determi- 
nation of topological charge. Spectral flows of the 
hermitian Wilson Dirac operator were also stud- 
ied to obtain some understanding of the role of 
topology in the mesonic spectrum in quenched 
QCD jn]]. Distribution of real eigenvalues of the 
massless Wilson-Dirac operator (identical to the 
zero eigenvalues of the Hermitian Wilson-Dirac 
operator) has received a lot of attention recently 
motivated by the need to understand the inability 
in extracting pion masses close to the chiral limit 
in quenched QCD 0. 
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